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Abstract 

We present a detailed analytical study of spherically symmetric self-similar solutions 
in the SU(2) sigma model coupled to gravity. Using a shooting argument we prove 
that there is a countable family of solutions which are analytic inside the past self- 
similarity horizon. In addition, we show that for sufficiently small values of the 
coupling constant these solutions possess a regular future self-similarity horizon and 
thus are examples of naked singularities. One of the solutions constructed here has 
been recently found as the critical solution at the threshold of black hole formation. 

1 Introduction 

In this paper we continue our investigations, started in |3J (referred to as I), of wave maps 
coupled to gravity, that is solutions of Einstein's equations with an SU(2) sigma field as 
matter. We found numerically in I that for a < 1/2 (a is the dimensionless coupling 
constant) the model admits a countable family of continuously self-similar (CSS) solutions, 
labeled by an integer nodal index n — 0,1, ... , that are analytic inside the past light cone of 
the singularity. We also provided evidence that the nth CSS solution can be continued up 
to the future light cone of the singularity if a < a n , where {a n } is an increasing sequence 
bounded above by 1/2. The purpose of this paper is to make the results of I into theorem- 
proof rigorous mathematics. This is accomplished by applying a shooting argument to the 
resulting dynamical system. We note that the case a = was previously analyzed in 0. 

The physical importance of the CSS solutions considered here was discussed in I, in 
particular we conjectured that in a certain parameter range (cto < a < a±) the n = 1 
solution is a critical solution at the threshold of black hole formation. This conjecture has 
been recently confirmed in numerical studies of the critical behaviour |3j and in the linear 
stability analysis 0]. As far as we know, this is the only case where the existence of a 
self-similar solution, which was numerically found as the critical solution in gravitational 
collapse, has been established rigorously. 
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2 Setup 



For the reader's convenience we repeat from I the basic setting for the problem. Let 
X : M — > N be a map from a spacetime (M,g ab ) into a Riemannian manifold (iV, Gab)- 
Wave maps coupled to gravity are defined as extrema of the action 



S= I (-ri^ + L W M)dv g (1) 



M 



V 167rG 



with the Lagrangian density 



fl ah A B 

Lwm = — ^-g a d a X d b X Gab- (2) 

Here G is Newton's constant and fl is the wave map coupling constant. The product 
a = 4:nGfl is dimensionless. The field equations derived from (1) are the wave map 
equation 

U g X A + Ti c (X)d a X B d b X c g ah = 0, (3) 

where T^ C (X) are the Christoffel symbols of the target metric Gab and D g is the d'Alembertian 
associated with the metric g ab , and the Einstein equations R ab — \g a bR = 8iiGT ab with the 
stress-energy tensor 

T ab = fl (d a X A d b X B - l -g ab (g cd d c X A d d X B )" Gab- (4) 

As a target manifold we take the three-sphere S 3 with the standard metric in polar coor- 
dinates X A = (F, 0, $) 

G AB dX A dX B = dF 2 + sin 2 F {dQ 2 + sin 2 9 d<$> 2 ). (5) 

For the domain manifold we assume spherical symmetry and use Schwarzschild coordinates 

g ab dx a dx h = -e- 2S A dt 2 + A~ l dr 2 + r 2 (d6 2 + sin 2 # d(j) 2 ), (6) 

where 5 and A are functions of (t, r). Next, we assume that the wave maps are corotational, 
that is 

F = F(t,r), & = e, $ = 0. (7) 
Equation (3) reduces then to the single semilinear wave equation 

n.F- -52Q = o, (8) 



where 



n g = -e 5 d t (e 5 A- l d t ) + ^d r (r 2 e- s A d r ), (9) 
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and the Einstein equations become 



d t A = -2arA(d t F)(d r F), (10) 
d r 5 = -ar((d r F) 2 + A- 2 e 2S (d t F) 2 ), (11) 

d r A = l -^-ar(A{d r F) 2 + A- l e 2 \d t F) 2 + 2^^\. (12) 
r \ r 2 J 

These equations are invariant under dilations (t, r) — > (At, Ar) so it is natural to look 
for continuously self-similar (CSS) solutions, that is solutions which are left invariant by 
the action of the homothetic Killing vector K — td t + rd r . To study such solutions it is 
convenient to use similarity variables p = r/(—t) and r = — ln(— t). Then K = —d T , so CSS 
solutions do not depend on r. Assuming this and using an auxiliary function Z = e s p/A, 
we reduce equations (8-12) to the system of ordinary differential equations (where prime 
is d/dp) 

F« + V-ad.Z^--^^^, (13) 

A' = —2apAF' 2 , (14) 
pZ' = Z(l + a(l-Z 2 )p 2 F' 2 ), (15) 

pA' = l-A-a(p 2 A(l + Z 2 )F' 2 + 2sm 2 Fy (16) 

The combination of (14) and (16) yields the constraint 

1 - A - 2asin 2 F + aAp 2 F' 2 (I - Z 2 ) = 0. (17) 

This system of equations has a fixed singularity at the center p = and moving singularities 
at points where Z(p) = ±1 and/or A(p) = 0. In terms of the similarity coordinate p, the 
metric (6) takes the form 

ds 2 = A-\l - Z' 2 ) p 2 dt 2 + 2A~Hp dtdp + A~H 2 dp 2 + t 2 p 2 [dQ 2 + sin 2 fl d<f) 2 ), (18) 

hence the hypersurfaces Z = ±1 are null (provided that A > 0). The first p\ where 
Z(pi) = 1 is the locus of the past light cone of the singularity at the origin (t — 0,r = 0) 
(in what follows we shall refer to the past and future light cones of the singularity as to 
the past and future self-similarity horizons (SSH)). By rescaling, p — > p/pi, one can always 
locate the past self-similarity horizon at p\ — 1, that is Z{1) — 1. To ensure regularity of 
solutions in the interval < p < 1, the equations (13-17) must be supplemented by the 
boundary conditions at both endpoints 

F(0) = 0, F'(0) = c, Z(0) = 0, A(0) = 1, (19) 
F(l) = |, F'(l)=b, Z(l) = l, A(l) = l-2a, (20) 

where c and b are free parameters. 

Our main result is the following theorem: 
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Theorem 1. For any < a < 1/2 and any nonegative integer n, the equations (13-17) 
have an analytic solution (F n ,A n ,Z n ) which satisfies the boundary conditions (19-20) and 
has precisely n oscillations of F n (p) around ir/2. 

In the next section we shall prove this theorem using a shooting technique. The numer- 
ical evidence for Theorem 1 was given in I . The case a = was proved previously in 
so hereafter we assume that < a < 1/2. 

3 Proof of Theorem 1 

For convenience we rewrite equations (13-15) in terms of H = F — tt/2: 

E" + ^-a(l + Z^ + A ^ H) z2) =0, (21) 

A' = -2apAH ,2 } (22) 
pZ' = Z(l +a(l - Z 2 )p 2 H' 2 ), (23) 

The constraint becomes 

1 - 2a - A + 2a sin 2 H + aAp 2 H' 2 {\ - Z 2 ) = 0. (24) 

The initial conditions at p = are 

H(0) = ~, H'(0)=c, A(0) = 1, Z(0) = 0, Z'(0) = 1. (25) 

Note that the above equations have a residual scaling symmetry p — > Xp. The initial 
condition Z'(0) = 1 is imposed temporarily in order to fix the scale. We shall refer to 
solutions of equations (21-24) satisfying the initial conditions (|25|) as to c-orbits. In the 
appendix we show that c-orbits exist locally and are analytic in p and c. Now we shall 
show that c-orbits can be extended up to a point p\ at which Z(pi) = 1. 

Proposition 2. For any < a < 1/2 and c > there is a pi(c) G (y/1 — 2a, 1), such that 
the c-orbit is defined for all p < p% and\im p ^ pi Z{p) = 1. Furthermore, the following limits 
exist 

-l < H= f lim H(p) <-, A = f lim A(p) = 1 - 2a cos 2 H, lim (1 - Z 2 )H' 2 = 0. 
2 p^pi 2 p^pi p^pi 

Proof: Let the maximum domain of definition of the c-orbit be < p < p\ and assume 
that Z(p) < 1 in this interval. Then, from constraint (|24] ) we have A > 1 — 2a > and 
hence A = lim p ^ Pl A(p) > (A exists since A(p) is monotone decreasing). By ( p3l) Z' > 0, 
hence Z = lim p ^ pi Z(p) exists. If Z < 1, then from constraint (^4j) H' 2 is bounded so 
H = lim p ^ pi H(p) exists, which in turn implies, again by (^4|), that lim p ^ pi H' exits. Thus, 
H, H', A, and Z all have finite limits at pi and therefore the c-orbit may be continued 
beyond p\ contradicting the maximality of p±. We conclude that Z = 1. 
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Now, we must show that H G (— 7r/2,7r/2) exists. Since Z — 1, we may no longer 
conclude that H' 2 is bounded but from equation (^) we get (In A)' = —2apH' 2 , so H' 2 is 
integrable near p\ which implies that H' is absolutely integrable (\H'\ < 1 + H' 2 ) and thus 
H exists. From constraint (|24[) , -ff(p) = ±7r/2 for some < p < p\ is not possible since 
1 - A > 0. Thus, -7r/2 < #(p) < tt/2 and so -tt/2 < H < tt/2. In fact, for p > p x /2 we 
have 1 — A > a > 0, so 2a cos 2 H > o > (remember that we assume a > 0), hence if is 
uniformly bounded away from ±7r/2, and thus —tt/2 < H < tt/2. 

To prove A = 1 — 2a cos 2 if, note that by ( p4[) <i = lim^^ H' 2 (l — Z 2 ) exits and is 
finite. Hence, by ( |2"BD lim p _ +pi Z' exists and is finite so 1 — Z 2 = 0(p — pi) near p\. If d ^ 0, 
then if'(p) ~ ^/(Pi — /°) would not be integrable near pi, thus <i must be zero. Inserting 
this into (H) we get A = 1 - 2a cos 2 H. 

Next, (Z/p)' > by ( p3|) and lim p ^ (^/p) = 1 by L'Hopital's rule, hence Z > p for all 
p > 0, and thus pi < 1. Finally, from ([^) and (|23D 

AZ 2 \' 2Z A AaH' 2 . . 

- < 0, (26) 



P 2 J P 

and since \im p ^ (AZ 2 / p 2 ) = 1, we have (AZ 2 /p 2 ) < 1 and hence pi > vA > \/l — 2a. 
If ^(p2) = 1 for some p2 < pi, we replace pi by p2 in the above arguments. 

Corollary 3. The function pi(c) is continuous. 

Proof: Let c be given and let e > 0. By Proposition 2, pi(c) is defined. The function Z(p) 
is monotone increasing for p < pi(c), so Z(px(c) — e, c) < 1, hence for all c sufficiently close 
to c, Z(px(c) — e, c) < 1, and thus pi(c) > pi(c) — e. To show that pi(c) < pi(c) + e for all 
c sufficiently close to c, we assume otherwise and get a contradiction. By the mean-value 
theorem Z(pi(c) + e, c) — Z(p, c) = Z'(£, c)(pi(c) + e — p). By continuity Z(p, c) is close to 
Z(p, c) and Z(p, c) is close to 1 if p is close to pi(c), hence Z(p, c) is arbitrarily close to 1. 
But Z'(p, c) > Z(p, c)/p > 1, so Z(pi{c) + e, c) > Z(p, c) + e > 1, which is a contradiction. 
Thus, pi(c) < pi(c) + e. 

Lemma 4. H'(p) is bounded near p\ if and only if H = 0. 
Proof: Suppose that H ^ and i?'(p) is bounded. Then, in (21) we have 

ttii i i i sin2iJ d 

H = bounded terms - - — - ~ , (27) 

Ap 2 (l-Z 2 ) pi-p 

where d ^ 0. This contradicts that H'(p) is bounded near pi and concludes the "only if 
part of Lemma 4. 

Suppose now that H(pi) = and H'(p) is unbounded. Without loss of generality we 
consider the case that H(p) < and H'(p) > near p\. Dividing equation (j2~H) by if' and 
integrating from p to pi we obtain 

pl / H" 2 , 9N , 2 sin(2iy) \ , 
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The first integral is divergent because lim p ^ pi In if' = oo. The second and the third terms 
are integrable (remember that H' 2 is integrable). Thus, to get a contradiction it suffices to 
show that the last term is integrable. We write this term as 

sin(2ff) = sin(2ff) H 
H'Ap 2 (l-Z 2 ) HAp 2 (l-Z 2 )H r { ' 

The first factor is continuous and we now show that the second factor is also continuous. 
Applying L'Hopital's rule we get 

H H' 1 

(1 - Z 2 )H' ~ p™ -2ZZ'H' + (1 - Z 2 )H" ~ p™ -2ZZ' + (1 - Z 2 )H"/H r ^ > 



Next, using (|2T|j we get 

~H" 2(1 -Z 2 ) , 9w 9x , 2 sin(2#) . , 

(1 ~ Z )jp= ~ + Ml + Z 2 )(l - Z 2 )H' 2 - -j^J-- (31) 

In the limit p — > pi, the first term on the rhs of ( |3~1~D obviously goes to zero, the second 
does by Proposition 2, and the third does by the assumption that H' — > oo. Thus, the 
limit (|30D is finite and consequently so is (p9|). This contradicts (p8|) and thus concludes 
the proof of the "if part of Lemma 4. 

Corollary 5. A c-orbit which has H(c)=0 is analytic on the whole interval < p < p\. 

Proof: The boundedness of H'(p) implies by ( plf ) that H" > —2H'/p is bounded below 
(remember that H(p) < and H'(p) > near pi), hence lim p _^ Pl H'(p) exists. Having 
that, it is easy to show by applying L'Hopital's rule to lim^^ii^^pi) — H^(p))/(pi — p) 
for k = 0, 1 that the solution (H, A, Z) is C 2 near p\. By a routine contraction mapping 
argument one can show that C 2 solutions are unique, hence a c-orbit must belong to the 
one-parameter family of analytic solutions from Proposition 14 (see the appendix). 

Next, we describe the behaviour of c-orbits for small and large values of the shooting 
parameter c. We define a nodal number of a c-orbit N(c) = number of zeros of the function 
H(p) on the interval < p < p\. We first show that c-orbits with small c have no nodes. 

Proposition 6. If c is sufficiently small then N(c) = 0. 

Proof: For c = we have H(p) = —n/2 and Z(p) = p so pi(c = 0) = 1. By continuity, 
for any e > and sufficiently small c we can find po such that 1 — e < po < Pi (c) < 1 and 
H{po) < — 7r/2 + e. We know from the proof of Proposition 2 that lim p ^ pi y/ pi — p H' = 0, 
hence 



f pl f Pl dp 

H(pi) — H(p ) = / H'(p)dp < const / —^=^= < const\fe. 

Jpo Jpo y/Pi-P 



(32) 

'Po -fpa VPl-P 

Thus, H(p) stays arbitrarily close to —n/2 all the way up to p\ if c is sufficiently small and 
therefore N(c) = 0. We remark that using a scaling argument one can derive the precise 
asymptotic behaviour of c-orbits for small c. We omit this argument since it is not needed 
for the proof. 

We show next that c-orbits with large c have arbitrarily many nodes. 
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Proposition 7. N(c) — > oo for c — ► oo. 

Proof: We rescale the variables, setting 

x = cp, #0) = H(p), A(x) = A(p), Z(x) = cZ(p). (33) 
Then, equations (21-24) become 

H" + -H>-a(l + ^)xH> 3 + „ M2H l =0, (34) 
* c Ar 2 (l-§) 

i' = -2asiiJ' 2 , (35) 
xZ' = Z(l + o(l - — )x 2 #' ), (36) 



with the constraint 



Z^ 

l-2a-A + 2a sin 2 # + aAx 2 H' 2 (l - — ) = 0, (37) 



c 2 



and the initial conditions at x = 



#(0) = ~, #'(0) = 1, 1(0) = 1, Z(0)=0, Z'(0) = 1. (38) 

As c — >• oo, the solutions of equations (34)-(38) tend uniformly on compact intervals to 
solutions of the limiting equations 

h" + h'-axh' 3 + S -^- = 0, (39) 
x ax z 

a' = -2axah' 2 , (40) 

xz' = z(l + ax 2 h' 2 ), (41) 

with the constraint 

1 - 2a - a + 2a sin 2 /i + aax 2 h' 2 = 0, (42) 
and the same initial conditions at x = 

h(0) = ~, h'(0) = 1, a(0) = 1, z(0) = 0, z'(0) = 1. (43) 



We observe first that the function a(x) is monotone decreasing by (40) and bounded below, 
a > 1 — 2a, by (|4~2"1) . Thus, no singularity can develop due to a going to zero. Also, by 
(fl2j) no singularity can develop due to /i' becoming unbounded. Thus, solutions exist for 
all x > (assuming the existence of a solution for small x). In order to complete the proof 
it is sufficient to show that the function h(x) has an infinite number of zeros for x > 0. 
Since a < 1, it follows from ( f42"|) that — 7r/2 < h(x) < n/2 for all x > 0. To show that /i(a;) 
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oscillates around zero we consider three cases: 

(i) Assume that linx^oo h(x) does not exist. Then, there must be a sequence . . . Xk < yk < 
Xk+i < Vk+i < ■ ■ ■ such that h has a local minimum at Xk and a local maximum at yk- By 
fl39f)) h\xk) = 0,h"(xk) > imply that sm(2h(xk)) < 0, hence h(xk) < 0. By a similar 
argument, /i(yfe) > 0. Thus, has a zero in each interval Xk < x < yk- 

(ii) Assume that a nonzero lim^oo h(x) exists. Then, from (f42|) lim^oo x 2 /i' 2 exists and, in 
fact, equals zero because lim^oo h(x) exists. This implies by ( |3TS| ) that lim^oo x 2 h"(x) = 
— sin(2/i(oo))/A(oo) ^ 0, hence lrnx^oo x 2 h' 2 (x) ^ 0. Thus the case (ii) does no arise. 

(iii) Assume that lim^oo h(x) = 0. We define the rotation function 9(x) by 

tanfl(x) = ^^, 0(0) =0. (44) 

Remark 1. The rotation function 9{x) is well defined because the case h(x) = h'(x) = 
is impossible for solutions satisfying the initial conditions ([43]) . To see this, assume that 
h(xo) = h'(xo) = for some xo > 0. Then, by (^) a(xo) = 1 — 2a and the unique 
solution with these initial conditions at x is h(x) = 0, h\x) = 0, a(x) = 1 — 2a for all x, 
contradicting the initial conditions ([4~3l). 

We want to show that lim^^oo 9(x) = — oo. Using (RPI) we obtain 



... s . 0n sin 2/i 2 cos 9 (l — 2acosh)sm9cos9 , 

x6>'(x) = - sir 6 1 . (45) 

2h a a 

Under the assumption lim^oo h(x) = 0, it follows from (|42|) that lim x _, OQ a(x) = 1 — 2a, 
hence for sufficiently large x 

1 / 2cos 2 0\ 3 
9'{x) w — sin 2 ^ + sin0cos^ + — < -— , (46) 

so lim^^oo 9(x) = — oo. Thus, given any integer there exists an Xk such that h(x) has at 
least zeroes for x < Xk- By continuous dependence on initial conditions, we may choose 
c > Xk/ \/l — 2a so that the c-solution has /c zeroes also for x < Xk- In terms of the variable 
p = x/c the c-solution has k zeroes for p < y/\ — 2a < pi(c). This completes the proof of 
Proposition 7. 

Next, we need two lemmas which tell us how the number of nodes N(c) changes under 
small variations of c. 

Lemma 8. If H(c) = 0, then N(c)=N(c) or N(c) = N(c) + 1 for c sufficiently close to c. 

Proof: First note that if H(p,c) has a zero at some p < pi(c), then H'(p ,c) ^ (see 
Remark 1) so by continuity of H(p, c) with respect to c, H(p, c) also has a zero if c is 
sufficiently close to c. Thus N(c) > N(c) and it suffices to show that N(c) < N(c) + 1. 
Let a < pi{c) be the last node of the c-orbit, that is H(a,c) = and, for concreteness, 
H (p, c) < for a < p < p\. By continuity with respect to c, H(p, c) will also have a zero 



S 



at a near a if c is near c. In order to prove that H(p, c) cannot have more than one zero in 
the interval a < p < pi(c), we now show that if H(p,c) becomes positive for some p > a, 
then it would not have time to change the sign again before going singular. Assume for 
contradiction that there is a segment a < pn < p < Pd of the c-orbit in which the function 
H(p) is monotone decreasing from a local maximum H(pn) > to H{pu) = 0. 
We define 

W = X -p 2 AH ,2 {\ - Z 2 ) + sin 2 H. (47) 



From @|) W = (A - 1 + 2a)/(2a), hence by © W < 0. We have 

H' 2 2 v 7 ^ , . 

W-sin 2 H pM(l-Z 2 )' - sin 2 tf p^A(l-Z 2 )' K ' 

Integrating the left-hand side from pn to pr>, we get (using Hn = H(pn)) 



PD -H'dp [ Hn dH [ Hn dH 



PN 
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> / > _ ) (49) 

y/W - sin 2 H Jo \jw - sin 2 H Jo \/sin 2 H N - sin 2 F 2 ' 

where the first inequality follows from W(p) < W(pn) = sin 2 Hn (since W decreases) 
and the second inequality is a simple calculation using a substitution sin if = u sin Hn 
(remember that Hn < tt/2). 

Next, we derive an upper bound for the integral of the right-hand side of ([481) . We have 

o° dp < 1 p dp < 1 p rfp 

We showed above that Z' > 1, hence 1 — Z > p x — p. Therefore 

/ 75 = < / = 2(Vpi - Piv - VPi - Pd) < 2Vpi - PiV- (51) 

J PN yl — Z J PN vPi-p 

By continuity of solutions with respect to c and by Corollary 3, pn is arbitrarily close 
to pi(c) if c is sufficiently close to c, hence it follows from (51) that the integral of the 
right-hand side of fl4"8|) is arbitrarily small. This is in contradiction with (|49|), hence H(p, c) 
cannot have a second additional zero, which completes the proof of Lemma 8. 

Lemma 9. If H(c) ^ 0, then N(c) = N(c) for c sufficiently close to c. 

Proof: Without loss of generality we assume that H(c) < 0. As above let a < pi(c) be the 
last node of the c-orbit, that is if (a, c) = and H(p,c) < for a < p < p\. Let a be the 
corresponding zero of H(p, c) for c near c. We want to show that H(p, c) cannot have an 
extra zero for p > a. Suppose for contradiction that H (b, c) = for some b > a. Then, 
there must be a 5 < b such that H (5, c) = H{c). Let us integrate the identity 

^ (52) 



y/W - sin 2 H Py/A(l - Z 2 ) 
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from 5 to b. For the left hand side we get 



H'dp r R dH 

(53) 



y/W - sin 2 H Jo y/\V -sin 2 H' 

From Proposition 2 we know that lim p ^ Pl (l - Z 2 )H' 2 = so W(p, c) < (1 + e/2) sin if 
for p near pi and hence W(p,c) < (1 + e) sin 2 if for c near c. Since is decreasing, 
W(S, c) <W(p,c) < (1 + e) sin 2 F. Thus 

B dH r R dH f 1 \ 7T . , 

> / ; == > arcsin —== > - (54) 



- sin 2 # Jo + e) sin 2 # - sin 2 // VVl + e/ 2 

for sufficiently small e, where the last but one inequality can be seen by substituting 
sin if = u sin H into the integral. By the same argument as in the proof of Lemma 8, the 
integral of the right hand side of (|52|) is 0(y/pi — p). By continuity of solutions with respect 
to c and by Corollary 3, 5 is arbitrarily close to pi(c) if c is sufficiently close to c, hence 
the integral of the left hand side of equation (52) is arbitrarily small. This contradicts (54) 
and completes the proof of Lemma 9. 

Now we are ready to make a shooting argument. We define a set 

C = {c\ N(c) = 0} (55) 

and let Co = sup Co. The set Co is nonempty (by Proposition 6) and bounded above (by 
Proposition 7) so Co exists. We claim that the co-orbit has no nodes and satisfies the 
boundary condition H(cq) = 0. To see this, note that the co-orbit cannot have a node 
because then by Lemmas 8 and 9 all nearby c-orbits would have a node so there would be 
an interval around Co without any elements of Co in it, contradicting the assumption that 
c is the least upper bound. Thus, iV(c ) = 0. Now, if H(c ) < 0, then by Lemma 9 all 
nearby c-orbits would have no nodes, so there would be an interval around Co consisting 
of elements of C , contradicting the assumption that c is an upper bound of C . Thus 
H(co) = 0. 

Next, we define C\ = {c > c | N(c) = 1}. This set is nonempty by the previous step 
and Lemma 8 and bounded above by Proposition 7, hence c\ = supCi exists. By the 
same argument as above, the Ci-orbit has exactly one node and satisfies H(c{) = 0. The 
construction of subsequent c n -orbits proceeds by induction. 

Conclusion of the proof of Theorem 1: 

Returning to the original variable F(p) and rescaling p — > p/pi(c n ) we get the solution 
of equations (13-17) which satisfies the boundary conditions (19) and (20) and has exactly 
n intersections with the line F = tt/2. By Corollary 5 this solution is analytic in the whole 
interval < p < 1. 
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4 Beyond the past self- similarity horizon 



In this section we consider the behaviour of the CSS solutions of Theorem 1 outside the 
past SSH, in particular we ask the question: do these solutions possess a regular future 
self-similarity horizon? Note that p = oo corresponds to the hypersurface (t = 0, r > 0) 
so in order to analyze the global behaviour of solutions (for t > 0) we need to go "beyond 
p = oo" . To this end we define, after I, a new coordinate x by 

We also define an auxiliary function w(x) = 1/Z(p). In these new variables, the past SSH 
where w = 1 is at x = 0, while the future SSH (if it exists) is at some xa > where 
w(x A ) = -1. 

In terms of x and w, the equations (21)-(23) become autonomous (where now prime is 
d/dx) 

H" - 2awH' 3 + ^^- = 0, (57) 
A(w z — 1) 

A' = -2aAwH' 2 , (58) 
w ' = -i + a {l-w 2 )H' 2 . (59) 

The constraint (24) becomes 

l — 2a-A + 2a sin 2 H + aAH' 2 (w 2 - 1) = 0. (60) 

From (20) the initial conditions at x = are 

H(x) ~ bx, w(x) ~ 1 - x, A(x) ~ 1 - 2a - 2a(l - 2a)b 2 x. (61) 

We know from Theorem 1 that for each a < 1/2 there is an infinite sequence {b n (a)} 
determining solutions which are regular inside the past SSH, that is for all x < (note 
that p = corresponds to x = — oo). In I we showed that for x > the solutions starting 
from the past SSH with the initial conditions (^l|) tend in finite "time" to w = — 1 if b is 
small, or to w = +1 if b is large. After I we shall refer to these two kinds of behaviour 
as to type A and type B solutions, respectively. Now we want to show that the solutions 
of Theorem 1 are of type A (and therefore possess the future SSH) provided that a is 
sufficiently small. Unfortunately, the shooting argument gives us insufficient information 
about the parameters b n so we cannot apply the above mentioned result of I to determine 
the character of solutions outside the past SSH. Instead, we shall make use of the obvious 
fact that for a = all solutions are of type A. 

Lemma 10. For sufficiently small a the c n -orbits of Theorem 1 (rescaled so that p\{c) = 1) 
have \b n \ uniformly bounded above for all n. 
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Proof: It was shown in || (see Lemma 4 in that reference) that for a = the solution to 
equations (57)-(61) for x < must exit the strip \H\ < ti/2 if \b\ is too large, say |6| > B. 
By continuous dependence, the same is true for sufficiently small a. But from Proposition 2 
the c-orbits must stay in the strip \H\ < 7r/2 for all x < 0. Thus, \b n \ < B for small a. 

Lemma 11. If a solution to equations (57)-(60) has w(xq) < and A(xq) > 2/3 for some 
Xq, then there is xa > %o such that \im. x ^ XA w(x) = —1, i.e., the solution is of type A. 

Proof: By (58) A is increasing for w < 0. Thus, using equation (59) and the constraint 
(60) we get for x > x 

. , ox ,2 1 — A — 2a cos 2 H 11 . . 

w ' = -l + a (l - w 2 )H' 2 = -1 + < -2 + - < --, (62) 

hence w must hit —1 for some finite xa > %o- 

Proposition 12. The c n (a)-orbits are of type A if a is sufficiently small. 

Proof: For a = and any b we have w(x) = 1—x and A(x) = 1; in particular A(3/2) = 1 > 
2/3 and it?(3/2) = — 1/2 < 0. By continuous dependence on initial conditions there exists 
a 5(b) such that if a < 5(b) and \b - b'\ < 5(b), then A(3/2,b') > 2/3 and w(3/2,b') < 0. 
This implies by Lemma 11 that the solutions corresponding to such values of a and b' are 
of type A. By a standard theorem of advanced calculus there is a 5' > (independent of b) 
such that the solutions with a < 5' and \b\ < B are of type A. By Lemma 10 any c n -orbit 
has \b\ < B, so for a < 5' the c n -orbits are of type A. 

By a similar argument as in the proof of Proposition 2 one can easily show that the 
type A solutions are generically only C° at the future SSH (for isolated values of a there 
are solutions that go smoothly through the future SSH). In I we showed that the leading 
order asymptotic behaviour at the future SSH is (using y = xa — x) 

w~-l+y, A~A -2aA C 2 y\n 2 (y), H ~ H - Cy \n(y), (63) 

where A Q = 1 — 2acos 2 if ) C = sm(2H )/2A , and H is a free parameter. Using this 
expansion one can check that the curvature is finite as y — ► which means that the type 
A solutions are examples of naked singularities. 
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Appendix (local existence theorems) 

In || (Proposition 1) Breitelohner, Forgacs, and Maison have derived the following result 
concerning the behaviour of solutions of a system of ordinary differential equations near a 
singular point (see also for a similar result): 

Theorem (BFM). Consider a system of first order differential equations for n + m functions 
u = (ui, ...,u n ) and v = (v h v m ) 

where constants Aj > and integers p>i,Vi > 1 and let C be an open subset of R n such that 
the functions / and g are analytic in the neighbourhood of t = 0, u = c, v = for all c G C. 



Then there exists an n-parameter family of solutions of the system (34) such that 



Ui{t) = Ci + 0(t"), Vi{t) = 0{t Ui ), (65) 

where Ui(t) and Vi(t) are defined for all c E C, \ t\ < to(c) and are analytic in t and c. 

We shall use this theorem to prove existence of local solutions of equations (21)-(23) 
near the singular points p = and p = 1. 

Proposition 13. TTie equations (21) -(2,3) admit a two-parameter family of local solutions 
near p = 

= -^ + cp + 0(p 3 ), (66) 

A(p) = l-acV + 0(p 4 ), (67) 
Z{p) = dp + 0(p 3 ), (68) 

which are analytic in c, d and p. 

Proof: Using the variables 

If + 7r/2 , 1-4 Z 

wi = , w 2 = H , w 3 = — —, w A = — (69) 

p p 2 p 

we rewrite the equations (21)-(23) as the first order system 

pw[ = —wi + w 2 , pw' 2 = 2wi - 2w 2 + p 2 hi, 

pw f 3 = -2w 3 + 2aw\ + p 2 h 2l pw f 4 = p 2 h 3 , (70) 

where the functions hi are analytic near p = 0. Next, substituting 

wi = tti - vi, w 2 = u\ + 2vi, 

UJ3 = v 2 + a(u\ — 2v\ — 8Mit>i), Wi = u 2 (71) 
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we put flTOD into the form (|6 



P u 'i = P 2 fu pu' 2 = p 2 f 2 , 

pv'i = -3ui + p 2 #i, p^2 = -2^2 + P 2 92, (72) 

where the functions fi,gi are analytic in an open neighbourhood of p = 0,ui = c,u 2 = 
d, v i = for any c and d. Thus, according to the BFM theorem there exists a two-parameter 
family of solutions such that 

Ul = c + 0(p 2 ) } u 2 = d + 0(p 2 ), (73) 

vi = 0(p 2 ), v 2 = 0(p 2 ), (74) 

which is equivalent to (66)- (68). 

Proposition 14. The equations (21)-(23) admit a one-parameter family of local solutions 
near p = 1 

H{p) = b (p-l) + 0((p-l) 2 ), (75) 
A{p) = l-2a-2a{l-2a)b 2 (p-l) + 0({p-l) 2 ), (76) 
Z(p) = p + 0((p-l) 2 ) (77) 

which are analytic in b and p. 
Proof: We define the variables 

u = H', Vl = - H', (78) 

P ~ 1 

(1 - 2a) - A . . , 2 Z -I 

v 2 = ± 2a 1 - 2a)H' 2 , v 3 = - 1. 79 

p- 1 p- 1 

Then, the equations (21)-(23) take the form (using t = p — 1) 

tu' = tf, tvl = - Vi + t 9i , (80) 

where the functions / and g\ are analytic in an open neighourhood of t = 0, u = b, = 
for any b > 0. Thus, according to the BFM theorem there exists a one-parameter family 
of solutions such that 

u(t) = b + 0(t), v l (t)=0(t) } (81) 
which is equivalent to (75)- (77). 
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